Warp Features in DBI Inflation 
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In Dirac-Born-Infeld inflation, changes in the sound speed that transiently break the slow roll 
approximation lead to features in the power spectrum. We develop and test the generalized slow 
roll approximation for calculating such effects and show that it can be extended to treat order unity 
features. As in slow-roll, model independent constraints on the potential of canonical inflation can 
be directly reinterpreted in the DBI context through this approximation. In particular, a sharp 
horizon scale step in the warped brane tension can explain oscillatory features in the WMAP7 CMB 
power spectrum as well as features in the potential. Differences appear only as a small suppression 
of power on horizon scales and larger. 



I. INTRODUCTION 

In Dirac-Born-Infeld (DBI) inflation [TJ [2], transient 
but rapid changes in the sound speed leave their imprint 
as features on the power spectrum. For string-motivated 
DBI examples, such features might arise from duality cas- 
cades which impart steps in the warped brane tension 
[3j H] . Annihilation of brancs during DBI inflation has 
also been shown to lead to particle production and to 
the imprint of features on the warp [5]. More generally, 
within the context of effective field theory [5] it has been 
shown that a sharp step in the sound speed leads to os- 
cillatory features in the power spectrum of fluctuations 

Power spectrum features from sudden changes in the 
warped brane tension of DBI inflation are closely related 
to those from sudden changes in the potential for canon- 
ical single field inflation. Measurements of the CMB 
temperature power spectrum from WMAP place obser- 
vational constraints on the latter. Recently, the gener- 
alized slow roll approach (GSR) [H [5] has been used to 
extract model-independent constraints from the WMAP 
data on features as sharp as 1/4 of an efold [10l HI] . 
Even sharper features lead to highly oscillatory power 
spectrum features which can evade these constraints due 
to projection effects. Indeed there is a special case where 
a sharp step in the potential on scales near the current 
horizon can fit the WMAP data better that a smooth 
model in the acoustic regime [12] . 

The GSR approach remains valid for single field infla- 
tion with non-canonical kinetic terms 13], including DBI 
inflation, with a suitable reinterpretation of the source of 
deviations from slow-roll |14) . In this Paper, we develop 
the GSR approach for DBI inflation and show how ob- 
servational constraints on potential features translate to 
constraints on warp features. 

In fjnj we briefly review the phenomenology of DBI in- 
flation and the exact computation of its power spectrum. 



In i III we develop and test the GSR approach in the DBI 
context and establish the correspondence between poten- 
tial features and warp features. In §IV[ we consider the 
special case of a sharp step in the warp analytically and 



show that it can explain the WMAP data as well as a 
sharp step in the potential. We discuss these results in 

9VI 



II. DBI INFLATION 

We consider DBI inflation to be a phenomenological 
model with the Lagrangian density 



£ 



1 - Vl - 2X/T(ct>) T(<i>)-V(<j>), 



(1) 



where the kinetic term 2X = — V^V^. In braneworld 
theories that motivate the DBI Lagrangian, <j> determines 
the position of the brane, T(<j)) gives the warped brane 
tension, and V(4>) is the interaction potential. 

As a consequence of the non-canonical kinetic struc- 
ture, field perturbations propagate at the sound speed 



,X) 



y/l-2X/T(<l>). 



(2) 



The inflaton energy density and pressure can be ex- 
pressed in terms of the sound speed as 



p&,x)= (^-i) T{ct>) + V{4>), 
p{<j>,X) = (l-c.)T(<j>)-V(<f>). 



(3) 



Note that for X/T C 1, c, = 1 and the Lagrangian, p 
and p take on their canonical forms. 

For the background equations of motion, we take the 
acceleration equation 



Hn 

IT 



2c s 



(4) 



where X — H 2 (f)\ 

Friedmann equation H 2 = p/3, and the field equation 



N /2, the Hubble parameter satisfies the 



(f>NN 



1 



+ 3c s 4> N - c\ 



H 2 



+ -(l-c s ) 2 (l + 2c s ) H , 



(5) 



Here and throughout the subscript N denotes d/d In a, 
the subscript <f> likewise d/d<f>, and we choose units where 
Mpi = (87rG) -1 / 2 = c = Ti = 1. If warp and potential 
features are absent near the initial conditions, initial val- 
ues for {(f), (fix, H} can be set on the slow-roll attractor 



c s V<t> 
3 H 2 ' 



V 

3' 



(6) 



where we assume that the Va term dominates over T§. 
Given that we choose to solve Eq. Q and ([5 1, we must 
ensure that the Friedmann equation is exact y satisfied 
on the initial condition [TS] . This can be achieved by first 
choosing the initial (f>(Ni), then taking 



H<f> 



N 



v_Y± 

3 V 




-1/2 



(7) 
Since 



and calculating p, H exactly through Eq. (|3| 
<f>N = H(/)n/H, we now have a self-consistent set of ini- 
tial conditions {0, 4>n, H} at A/j. This technique remains 
valid for all c s in the slow-roll approximation. On the 
other hand, for c s <C 1 the slow-roll approximation can 
remain valid even for steep potentials in Eq. ([6]). 

We evolve these equations until the field reaches <fr = 
4> cn d which we take to be the end of inflation and define 
N = to be this epoch 



N = 



In a e 



In a 



din a, 



(8) 



such that N < during inflation. For the purposes of 
calculating the power spectrum, it is useful to express 
the efolding number, N, in terms of the sound horizon, 
the comoving distance sound can travel from N to the 
end of inflation 



s(N) = I dN ° 



1 



iV 



aH 



dN- C> 



iV 



H 



(9) 



By defining the effective reheat temperature as T re heat 
7o/a onc i, where the present CMB temperature is To 
2.726K, we can express the sound horizon as 



-65.08 



TV, 



-heal 



500Mpc JN 
The curvature power spectrum is then given by 



dN ak- 



(10) 



with ' = d/dhis throughout and 



f = 8ir 



2 t H c s ( aHsY 

"a J 



H 2 



(14) 



We solve Eq. ( 12 ) assuming Bunch-Davies initial condi- 
tions where lim^ s> i y = e lks . Note that written in this 
form, the Mukhanov-Sasaki equation takes exactly the 
same form for canonical and non-canonical kinetic terms. 
This fact will allow us to remap existing constraints on 
S'(lns) onto DBI models. 



III. GSR APPROXIMATION 

In this section, we begin by reviewing the GSR ap- 
proach to calculating the curvature power spectrum 
[H [HI Q3| and show how to apply it to the DBI model. 
We then test the accuracy of the approach against the 
exact computation for a step-like feature in the warped 
brane tension [4] and show that the model-independent 
constraints on features in the potential for canonical sin- 
gle field inflation [TU1 E] can be readily reinterpreted in 
the DBI context. 



A. Technique 

Briefly, the GSR approach to solving the Mukhanov- 
Sasaki equation ( 12 1 is to consider the RHS as an external 



source with an iterative correction to the field value y. 
To lowest order, we replace y — > yo where 



2/o= l 



kt 



(15) 



is the solution to equations with g — > and solve for 
the field fluctuation y through the Green function tech- 
nique. To second order in slow-roll, the curvature power 
spectrum is given by |10[ I17j 



lnA|, = G(ln.s min ) + J 



ds 



W{ks)G'{\ns) (16) 



A 2 _ k 3 Pn v 
A K = = lim 

27T 2 ks-+0 



ksy 



f 



(11) 



where the mode function y obeys the Mukhanov-Sasaki 
equation [14l H6] 



d*y 

ds 2 



k 2 - 



y = 



g(lns) 



Here 



9 = 



f" ~ 3/' 

/ : 



(12) 



(13) 



where the window function 

3sin(2u) 3cos(2it) 3sin(2u) 



W(u) 



Here 



and thus 



2u 3 u 2 2u 



G = -21n/+-(ln/)', 



G' = -2(ln/)' + -(ln/)" 



l^-flCm/) 



/l2 



(17) 



(18) 



(19) 
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We call G' the GSR source function. The quadratic term 
in (In/)' appears to ensure constant curvature fluctua- 
tions above the sound horizon. 

The I\ and / 2 integrals are the second order corrections 

h{k) = 4= / —G'(\ns)X{ks), 
V2 7 s 

h(k) = -A -[X + -X'] J -F 2 {u) 1 (20) 

Jo w 3 / 



with u = ks, 



and 



F 2 (w) 



« 2 /' 



X(u) = -^(sinu — itcosit)^ 



(21) 



(22) 



To calculate the power spectrum in the GSR approx- 
imation, we need to obtain the source functions G' and 
/'// from the solution to the background equations of 
motion Q and ([5|. In terms of the slow-roll parameters 



G' = 



2,„ „ x 2,aHs 1X , 

- (2e H - 2 m - a x + l) 2 

3 3 c s 

2 . ai7 s 

+-( 1)(4 + 2e H -2r, H -ai) 

3 c. 



(23) 



1 / flgs 

3 1^7 



2<5 2 + 2e|f - 2tj h - 2t& 



-3cti + 2r) H <r 1 + of - e H (Arj H + 01) - a 2 



and 
f 



aHs 
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Vh - £h+ 2 CT i 



1 - 



aHi 



(24) 



where the additional slow-roll parameters are defined by 

1 d In e^f 



»7h 


" £ff 2 


dN 
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= tHVH + 




dN 




din c s 
~ dN ' 






cr 2 


do - ! 
= dN' 







Using the field equation ([5]), we can write 



'WW 



■> N 



£H - c 2 s fj H 



1 (l-c.)(l + 2c s )_ 

2 



(25) 



(26) 



where 



m = 3 



(27) 



These auxiliary parameters t)h and eri quantify slow-roll 
deviations generated by features in the potential and 
features in the warp respectively. 

In terms of the auxiliary parameters, the slow roll pa- 
rameters themselves become 



1 



1 - c 



1]H 



-m 



2 ' 2 1 + c. 

„2\ ~ 



8 ~ 



(7i = (1 - c B )ai + (1 - c s )fy ff . 



(28) 



Note that for r/n, the term involving d\ is suppressed 
both as c s —¥ and c s — > 1. Furthermore rys is slow roll 
suppressed on the attractor of Eq. (|6| and for c s = 1, 
= Vh- If features in T(</>) drive deviations from slow- 
roll then the &\ term dominates, vh — [c s /2(l + c s )]ai 
and hence \rj H \ <§; |cri| for c s <C 1. 

The remaining slow roll parameters 02 and S 2 can be 
constructed by taking the derivatives of <j\ and vh 



0-2 = (1 - C s ) 

s 2 - 



da 1 
dN ~ 
l + c 2 s dfj H 



- 2c s Uif] H , 



£hVh + Vh~ c 2 s air) H 



2 dN 

c a 1 — c s dai c s 1 — 2c 



2 1 



2 (l + c s ) 



c s - 

3— ClCTl, 



(29) 



where 



4>N 



2e H 



d<j\ 



dN 



T 

T 4>4> 
T 



4>N 



T 



(30) 



(l + c. 5 + 2c 2 ). 2 



2 XX 



2(1 + c s ) 



<*\ + (£H - c s r] H )a 1 



For sharp features in the warp and potential, er 2 and <5 2 
dominate respectively in G' due to the appearance of sec- 



ond derivatives in Eq. (|30|). Note that since 4> 2 N = 2e#c s , 



the impact of fractionaTfeatures in the warp vs. the po- 
tential potential is suppressed by slow roll parameters. 

In the slow roll limit, one can iteratively substitute the 
attractor solution Eq. ([6| into the field equation to obtain 



G' w 4e H - 2r ]H + o x 

'V 2 
v 



c s (2 + c 2 ) 



2c 



y<f><p 

V 



Cs(l-C 2 s ) 



(31) 
T V ' 



where c s ((f>) is given by the attractor solution Eq. (J7|. 
The absence of a term in Eq. (31 ) can be attributed 



to the fact that the attractor solution is determined by 
Vf/,. Furthermore, in the slow roll limit, evolution in G' 
is second order in slow roll parameters and 



/' „ 1 

— ss -2<l h + Vh - -01 



1 



G' 



(32) 



3 



and so 



h 

h 



2V2 



-4 



(Ae H - 2ri H + 01), 
/ 



(4 eiT - 2r? H + aiY 



(33) 



Thus in the slow roll limit the total second order correc- 
tion involves a near cancellation of the I\ and I2 terms 



1 



1 



1 



In [1 + -If(k) + -I 2 (k)Y + -If(k) }*lf + h 



1 (4e H - 2 m + (TxY 



(34) 



If sharp features in either the potential or the warp 



dominate 



G' 



-a 2 



-2c, 



V 4>4> 
V 



1 — Cg 
3(1 + c s )^T 



b 2 N . 



(35) 



Unlike the canonical case, the functional constraints im- 
posed by observational constraints on V{4>) and T((f>) dif- 
fer in the two limits. Nonetheless deviations in the power 
spectrum for the sharp feature case, which can be large 
yet observationally viable, share strong similarities be- 
tween those generated by T((j>) and V(4>). Also, unlike 
the slow-roll limit, Ii tends to be larger than I 2 in that 
only it depends directly on second derivatives of T(</>) or 
V((f>) as we shall see explicitly in the next section. 



B. Numerical Tests 

The GSR construction in the previous section applies 
to any model with features in the warped brane tension 
T((p) or potential V((f>). For definiteness and motivated 
by the WMAP data, we test the GSR approximation on 
models where T(cf>) has a step feature [I] 



T(<P) 



-[1 + bF^)} 



with 



F((j>) = tanh 



- 1 



and <p inflates on a potential 



V(ct>) = V (l - ~P<P 2 ), 



(36) 



(37) 



(38) 



rolling from small to large values. We have chosen a con- 
vention that after the feature, T(</>) goes back to its b = 
value. Since physical scales are matched to the end of in- 
flation through Eq. ^ , this simplifies the comparison to 



the smooth featureless case. For simplicity, we will take 
T reW = Vq /4 following [18]. Thus the DBI step model 
is specified by 4 parameters {As, Vo, /?, </> e nd} controlling 
the underlying smooth spectrum and 3 parameters de- 
scribing the step feature {4> s ,b,d}. 

In order to set the parameters for the smooth b = 
spectrum, it is useful to re-express the attractor solution 
of Eq. ^ in terms of efolds [3] 



'■>N 



dN HJY, 



(39) 



and hence 



<P(N) 
Cs(N) 
e H (N) 



'N-N ' 



1 



P N - N : 



(40) 



where N is an integration constant determined by our 
definition that A(</> cn d) = 0. In the slow roll approxima- 
tion 



A 2 



( H 



H 2 



(N - Nq)* 
4tt 2 A b ' 



8ir 2 e H c s 



and the tilt 



din A 

rfln/c 



■R 



4 

N - N ' 



(41) 



(42) 



Note that to have a tilt that is compatible with ob- 
servations n s — 1 ~ —0.04 at N ~ —50 one requires 
N ~ 50. In string-inspired models where </> en d ~ Hy/Xs, 
No = 0(1). While such problems can be ameliorated by 
introducing stringy physics not included in the DBI ac- 
tion [18] doing so degrades the predictive power of cal- 
culations based on this action (c/. [3]). We therefore in- 
stead require e nd <S Hy/Xs so that inflation ends while 
the field is deep in the DBI regime. Finally, to satisfy 
constraints from upper limits on equilateral type bispcc- 
tra, we require c s > 1/30 for scales relevant to the CMB 

These conditions are satisfied by the following choices 

A B = 1.93 x 10 15 , 
V = 7.10 x 10~ 26 , 
P = 0.5, 



6 end = 1.065 x 10" 



(43) 



In our parameterization Vq drops out of expressions for 
the curvature power spectrum at a fixed N — No or n s — 1 
and only impacts the mapping between field and physical 
scale through T 1G hoat ■ It also enters into the tensor-scalar 
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FIG. 1. GSR source function G' for a warp step with b — —0.4 
(top) and b = -0.005 (bottom) with d = 2.81 x 10 -11 Also 
shown is the approximate form based on in the sharp, 
small amplitude limit from Eq. ( 45 1 which is an excellent ap- 



proximation for the small b case and remains in good quali- 
tative agreement for the high b case. 



ratio and so we take for definiteness a small value such 
that tensors are negligible. 

Since the GSR approximation reproduces the exact 
second order expansion in slow roll parameters by con- 
struction when they are all small, we test the technique 
for the nontrivial case where b is order unity. 

In Fig. [lj we show an example of the GSR source func- 
tion G' where b = —0.4 (top) and —0.005 (bottom) and 
d = 2.81 x 10" 11 and <j> s = 5.67 x 10~ 8 . Both cases ap- 
pear like the second derivative of the step in T(<j>) with 
a width determined by the number of efolds it takes for 
the inflaton to cross the step 



S In 



5N 



d 



(44) 



The main difference is that at the larger b value the lo- 
cation, amplitude and width of the feature differ slightly. 

In Fig. [2] we show the corresponding power spectrum. 
In the top panel, we compare the power spectrum from 
the full GSR approximation ("GSR2") of Eq. @ to the 
exact solution. In the middle panel, we show that the 
approximation is accurate at the 1-2% level for the or- 
der unity feature. Moreover, the second order correc- 
tions remain small as shown in the bottom panel where 
"GSR1" denotes setting 7 2 = and "GSR0" denotes set- 
ting both Ii — I 2 — in Eq. (16). Here the maximum 
value that |7i| attains is 0.37. As in the canonical case, 
max|7i| < l/v2 ensures accuracy in the power spectrum 
of the GSR approximation, typically to a few percent in 



10- 



0.0001 
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0.001 



0.01 



FIG. 2. GSR vs exact solution for the power spectrum 
(top panel), the fractional difference between the two (middle 
panel), and the impact of second order corrections (bottom 
panel) corresponding to the b = — OA model in Fig. [IJ (top) . 
GSR2 denotes the full solution {16} with 7 a and 7 2 , GSR1 the 
solution with 7 2 = and GSR0 with 7i = J 2 = 0. While 
the GSR2 solution captures effects at the 1 — 2% level for 
b = —0.4, even the leading order GSR0 is accurate at the 
10% level. 



observables such as the CMB power spectrum [TT] . 

Note 7 2 provides a negligible absolute correction for or- 
der unity and smaller features. For small features 7i and 
I 2 corrections do become comparable but in that case 
both are negligible (see Eq. 33 1. Since both leading or- 



der and 7i terms depend only on a single source function 
G'(lns), observational constraints from the power spec- 
trum may be directly mapped onto constraints on this 
GSR source function [3]. 



C. Observational Constraints on Broad Features 

Model-independent analysis of features in the source 
function G' have been conducted using a principal com- 
ponent basis that is complete for the WMAP7 data set 
for broad inflaton features that are traversed in 6 N > 1/4 
[TUJ [TT] . In the acoustic regime of s ~ 100—400 Mpc, con- 
straints on possible deviations are extremely tight with 
percent level limits on the broadest features [10]. As il- 
lustrated in Fig. [TJ these constraints can be interpreted 
in the DBI context as limits on 



G' 



2(1 -c.) 



3(1 



T <P<P 
;C s e H — 

c s ) T 



(45) 



The only broadband feature that marginally improves the 
likelihood is associated with the known I ~ 20 — 40 glitch 



in the WMAP data. However, a simple step in the warp 
does not fit the data as well as a step in the potential due 
to the change in c s across the step [I] . Since the attractor 



and so 



solution for the power spectrum in Eq. ( 41 ) depends on 
the sound speed, a step leads not only to oscillations, but 
also a step in the power spectrum across the feature (see 
Fig. [2}. 



IV. SHARP STEP 

Sharp features in the warped brane tension that are 
traversed by the inflaton in 8N <C 1 produce high fre- 
quency oscillations in the curvature power spectrum. 
These are more difficult to constrain observationally than 
broadband features due to projection effects and sky cov- 
erage. They are also significantly more cumbersome to 
calculate as their effects persist over orders of magnitude 
in wavenumber. 

In this section, we will derive an analytical solution for 
a very small and very sharp step in the warp factor and 
then test it against order unity steps. From Fig. [3J we 
can infer that for 6 <C 1 an analytic model based on inte- 
grating derivatives of T{cf)) should be accurate once the 
appropriate conversions between the field and sound hori- 
zon are made. At larger 6 we can see that the main differ- 
ences are that the location, amplitude, and width of the 
feature in G changes which we shall see require a recal- 
ibration of corresponding effects in the power spectrum. 
We use this analytic approximation to show that there is 
a DBI equivalent to the sharp potential step model that 
improves the WMAP7 likelihood by 2Aln£ ~ 12 [15] . 

We first start with some general considerations dic- 
tated by energy conservation and the slow-roll attractor. 
If the inflaton crosses a step in 8N <C 1 then we can 
ignore energy loss to the expansion and set the total en- 
ergy p in Eq. O) to be equal before and immediately after 
the crossing 4|. Kinetic energy in excess (or deficit) of 
the attractor after the step will then dilute away on the 
8N ~ 1 timescale. Denoting with A the change in quan- 
tities going through the step, we have immediately after 
the step 

Ac s = 1 - c s AT 
c s l + c s AT/T T ' 1 ' 

Note that for a decrease in T, energy conservation re- 
stricts an amplitude of \AT/T\ = \2b\ < l/c s . For a 
small amplitude warp feature, we can linearize 



(l-c s ) 



AT 
~T~' 



(47) 



Thus for the case of a small, sharp step in T, the sound 
speed takes a fractional step of comparable amplitude. 
Furthermore the slow-roll parameters a± and a 2 follow 
by taking derivatives of AT /T during the interval around 
the step. Similarly 
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P 



1 



T 

V' 



(48) 



Ae H 



1 



AT 



1 + c s AT/T 1 + c s 
c. Ac., 



T 



1 



(49) 



Note that at low sound speed, the relative effect of the 
step on tn is suppressed vs c s by c s /(l + c s ), as are t]h 
compared with a±, and 82 compared with 02, in agree- 
ment with Eq. 



(281 



After crossing the step, we know that the inflaton hits 
the attractor solution ([6| as the kinetic energy from the 
step decays after several efolds. For a small amplitude 
step 



Ac 



1 , „ AT 
c7 = 2 (1 - C ^' 



and 



Aeg 



Ac s 

Cs 



(50) 



(51) 



Given that the change in e# is determined by the change 
in the sound speed we seek to quantify the full evolution 
of c s from the step through to the attractor regime. 
Following [T23 , we begin by expanding the field as 

= ^> + 0i, (52) 

where tf> = <f>o when 6 = 0, and calculate to zeroth order 
in the unperturbed slow-roll parameters. Furthermore, 
the expansion rate is also unaffected by the warp feature 
since ejj <C 1 throughout. The field equation for <f>\ then 
becomes 
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2 u T i 

c sO)<PN0 7f- 



l (l-c s0 ) 2 (l + 2c a0 ) 
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"sO 



T41 
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(53) 



where again and 1 denote unperturbed 6 = and finite 
6 perturbations respectively. Here we have used the fact 
that 



Csl 
CsO 



zc s0 



■>N1 



Ti 
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(54) 



We can further transform the time variable from efolds 
N to background field value <po by taking <fipfo ps const. 
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(55) 



The first term on the RHS can be integrated by parts 
to make the whole source proportional to Tim. For sharp 
features, Ti^ is very concentrated around the feature and, 
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FIG. 3. Analytic vs. GSRO solution for a small amplitude sharp step b = —0.005, d = O.OO50jvo = 2.44 x 10~ 12 , with 
c„o = 0.0507 (left) and 0.50 (right). Top panel: difference in In Pn between this model and the same 6 = model. Bottom 
panel: difference between the curves in the top panel divided by the smooth envelope of the oscillations (see text). 



consequently, we can approximate the background quan- 
tities by their values at 4> s . Combined with the boundary 
condition that the field is on the attractor before the step 



0<0 S dcj)() 



(56) 



we obtain 



d<t>i 
d<f> 



1 



-s0 



C s0 1 ~ c s0 

2 l + c s0 



3(4. s -4> ) 



b[F{(j> ) + 2}e *«o 



Using this result in Eq. (54) and replacing 



(f> s -<f> = (f>No(Ns-N), 



(57) 
(58) 



we obtain 

c,i _ 1 - c 2 

CsO 



(59) 



Note that before the step c s i/c s o = — (1 — c 2 )6 and right 
after the step c s \/c s q = (1 — c s o) 2 b and so Ac s /c s — 
2(1 — c s0 )b as expected from Eq. (47). Several efolds 
after the step c s \/c s q = and so Ac7/c s = (1 — c 2 )6 as 
expected from Eq. (]50h. 



(60) 



2 (l + c s0 ) 



From en = <j> N /2c s we obtain 



em 



-WW 



which also satisfies expectations from Eqs. (49 1 and (51 1 . 

From these quantities, we calculate G taking e# — > 0, 
b -> 0. d -t 



G' « --a-a + -<5 2 - -tri - 2r? H + 1). (61) 

3 3 3 3 c s 



In this approximation, 



aHs 



-Is- 



5 - 2c 



s() 



3c : 



.sO 



N—N* 



bF(cf> )e 

v s o 

- 3(1 ~ Cs0) V (^o) + ^e 3 ^-"^. (62) 
8 

After several integ rati ons by parts we obtain the change 
in In A 2 ^ from Eq. ( 16 1 due to the feature, with 1\= 1-2 = 
from the feature, as 

In A 2 ^ = C 1 W(ks s ) + C 2 W'{ks s ) + C 3 Y(ks s ), (63) 

where 

. . 6a;cos(2a;) + (4x 2 - 3)sin(2x) 
Y(x) = - 3 (64) 

is proportional to J dlnxW' /x. Here 

d = 2(1 - c 2 )6, 
21-c s0l 



C 2 

c 3 



3 1 + c s0 
5 — 2c s0 — 3c s0 
4 



(65) 
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FIG. 4. Analytic vs. GSRO solution for a large amplitude sharp step 6 = —0.25 and d = 2.44 x 10 12 with low sound speed 
CsO = 0.0507 model as in Fig. [3] Left: the analytic approximation using the linear in b scalings of Eq. (651 begins to depart 



from GSRO as b approaches unity. Note however that the functional form of the feature remains the same but amplitudes 
and damping require rescaling. Right: nonlinear ly rescaled amplitudes and damping of Eq. \72\ recover the few percent level 
accuracy seen in the small step case. 



Given that 



lim W[x) 

x — ^0 

lim W'(x) 

X— ^0 

lim Y(x) 



1, lim W(x) = 0, 

x—>oo 

0, lim W'ix) = -3cos(2a;), 

0, lim Y(x) = 0, (66) 



the W term represents a step in the power spectrum at 
ks s ~ 1 of fractional amplitude 26(1 — c^q), which follows 
directly from the attractor solution, and the W term 
represents a constant amplitude oscillation at ks s 3> 1. 
The latter is exactly the same form as oscillations pro- 
duced by a step in the potential for a canonical kinetic 
term (see [12] Eq. 32). Unlike the canonical case, the step 
in power is comparable to the amplitude of oscillations. 
Furthermore, the additional Y term changes the solution 
near ks s ~ 1. Since even for fc< 1 a small error in the 
location of the feature s s , which controls the frequency 
of the oscillation, causes a noticeable change in the phase 
of the oscillation over many cycles, we define it such that 



8G'Qns a ) = 



(67) 



for the change from the smooth 6 = model. This def- 
inition differs slightly from the sound horizon at <fi s for 
large 6 as shown in Fig. [l] 

For finite step width d in field space, the infiaton tra- 
verses the step in As/s s ps |dln s / d<fi\d w d/^jvo- The 
window functions W and W' oscillate on a time scale 



As = 1/k. Thus the integral over G' is damped for 
ks s > 4>N/d. For the tanh step, the integral can be 
approximated following |12j 



lnA^ 



CxWikss) + C 2 W'(ks s ) + C 3 Y(ks s ) 
ks s 

Xd 



x V 



where 



1 



Xd 



d\ns nd 
and the damping function is 

T>(y) 



?N0 

nd 



V2i 



H0 c s0 



nd 



V 

sinhy 



(68) 



(69) 



(70) 



To obtain the full power spectrum we add lnA^ to a 
calculation of the 6 = model. This can be an exact nu- 
merical solution, a slow-roll approximation, or the GSR 
approximation. For comparison purposes, we choose here 
to take the GSRO (h 



III 



I2 = 0) solution from 
In Fig.[3j we test the analytic approximation for a small 
amplitude sharp step 6 = —0.005 and two values of the 
sound speed. In the lower panel we divide the difference 
between the analytic and GSRO solutions by the envelope 
function 



3C 2 V 



ks s 

Xd 



(71) 
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/ 



FIG. 5. CMB temperature power spectrum for the a horizon 
scale, sharp step in the potential V(<f>) of a canonical inflation 
fit to the WMAP7 data versus the matching step in the DBI 
warp T((f>). Lower panel shows the fractional deviation from 
a smooth spectrum with the same average power. The two 
types of steps produce nearly identical deviations at high £ 
and so fit the data equally well. Both spectra are calculating 
using their respective analytic approximations and without 
gravitational lensing. 



tor solution and the form of the C2, C3 corrections is 
motivated by the fact that T = before the feature for 
b = 1/2. In Fig. E] (right) , we show that the agreement is 
again good after these rescalings even for the b = —0.25 
case. 

Given that the analytic approximation works quite well 
even for relatively large values of b and its functional form 
mimics a step in the potential at ks s 3> 1, we can remap 
results for the latter onto the former. A potential step 
at s s — 8163 Mpc of amplitude C2 = 0.11 improves the 
WMAP7 likelihood by 2Aln£ « 12 QJ]. In terms of the 
warp step, these parameters translate into b = —0.218 
for c s q = 0.0507 and d 0. In Fig. [5] we compare the 
CMB temperature power spectra predicted by the two 
models using the best fit parameters for the cosmological 
parameters: fl b h 2 = 0.0222, fl c h 2 = 0.11, h = 0.71, 
t = 0.10, lnl0 10 A s = 3.077, n s = 0.965 such that the 
underlying smooth power spectrum is 

a ™=Mos^r" (73> 

Note that aside from small changes at low multiple £ 
where the cosmic variance is high, the two spectra are in- 
distinguishable. Thus a step in the warp fits the WMAP7 
data as well as a step in the potential. 



Agreement is at the 1% level except on scales much larger 
than the step ks s <C 1 and those affected by damp- 
ing ks s > Xd- In the former case differences from the 
change in (f>N from (/>tvo due to the different slow-roll at- 
tractor change the mapping between <fi and Ins. Near 
the damping scale, small changes in x<i are amplified in 
the fractional difference due to the exponential nature of 
the damping even though the absolute prediction remains 
accurate. 

As b approaches order unity, the analytic approxima- 
tion begins to misestimate the amplitudes and damping 
of the features. In Fig. [4] (left) we show a case where 
b = —0.25 and d = 2.44 x 10~ 12 , with low sound speed 
c s q = 0.0507. Note that the functional form of the power 



spectrum in Eq. ( 63 1 remains the same only the coeffi- 
cients differ. We therefore rescale them as 



Ci = In 

C 2 ■■ 



1 - 26c^ 
1 - 2b 
2 1-c.o 1 
3 1 + c s0 Vl - 26 



C, = 



Xd 



5 - 2c 



s0 



3c 



s0 



1 



^/l - 26 



b, 



do 



d In ; 



1 

ird 



(72) 



where the form of C\ can be derived from the attrac- 



V. DISCUSSION 

We have shown that the GSR approximation can be 
applied to DBI inflation to constrain features in the 
warped brane tension T(<p) from observational data. The 
approximation accurately recovers corresponding fea- 
tures for up to order unity deviations. Previous work on 
constraining the GSR source function G' and hence sec- 
ond derivatives of the potential V(<j>) for canonical fields 
can be directly reinterpreted in the DBI context as lim- 
its on the second derivative of T(<j>) [Tffl 111) . The main 
difference between the two is that features in T(</>) once 
traversed can strongly affect the slow-roll attractor for 
modes that cross the sound horizon later. 

The correspondence between features in V(4>) and 
T(cj)) is especially close in the limit of extremely sharp 
features, for example a step feature. In both cases the 
power spectrum exhibits constant amplitude oscillations 
for modes that cross the sound horizon after the step. 
Consequently, the preference for a horizon-sized step in 
the potential in WMAP7 implies a corresponding prefer- 
ence for a step in T(0). The main difference is a reduction 
of power for the low k modes that cross before the fea- 
ture. The large cosmic variance of these modes prevents 
a significant distinction between the two. On the other 
hand, features in V(<p) for canonical inflation and T((j>) 
for DBI inflation should induce very different bispectra. 
We leave these considerations to a future work. 
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